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Abstract The rotating magnetic B stars are a class of vari- 
ables consisting of He-strong and some P Cep stars which have 
oblique dipolar magnetic fields. Such stars develop co-rotating, 
torus-shaped clouds by channeling wind particles from their 
magnetic poles to circumstellar regions centered around the 
plane of their magnetic equators. The rotation of the cloud-star 
complex permits the study of absorptions from the cloud as it 
occults the star. In this paper we describe a quantitative analysis 
of archival lUE data to map the properties of these clouds over 
four stars (HD 1 84927, a Ori E, P Cep, and HR 6684). By com- 
puting spectral synthesis models for these stars, we find that 
only P Cep has a solar-like metallicity. Our analysis also shows 
that the metal composition across the surfaces of all these stars 
is at least approximately homogeneous. 

Using the Hubeny code CIRCUS, we demonstrate that the 
periodic variations of broad-band ultraviolet continuum fluxes 
can be explained fully by the absorptions of the co-rotating 
cloud. We show next that among selected lines, those arising 
from low-excitation states are selectively affected by cloud ab- 
sorption and turbulence. Our analysis also quantifies the cloud 
temperatures and column densities required to match the ab- 
sorptions of a number of weak to moderate strength resonance 
lines. These temperatures increase with the ionization poten- 
tial of the parent ions of these various lines, a result which is 
consistent with radiative equilibrium models in which temper- 
ature increases with proximity to the star's surface. Although 
these attributes appear stable from one epoch to another, dy- 
namic processes are nonetheless at work. Both the strengths 
and widths of resonance lines at occultation phases indicate the 
presence of a turbulence in the cloud which increases inwards. 

The spectroscopic hallmark of this stellar class is the pres- 
ence of strong C IV and N V resonance line absorptions at oc- 
cultation phases and of redshifted emissions of these lines at 
magnetic pole-on phases. The emissions have characteristics 
which seem most compatible with their generation by high- 
energy shocks at the wind-cloud interface, as predicted recently 



by Babel (1998) 



Key words: stars: individual pCephei - stars: individ- 
ual HR6684 - stars: individual oOriE - stars: individual 
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1. Introduction 

The He-strong Bp stars are a subclass of early-type B stars 
with an array of remarkable properties, chief among which are 
their oblique dipolar magnetic fields, magnetically controlled 
winds, and chem ical s urface anomalies (e. g. Boh lender et al. 



1987, Bohlender 1994 Bolton 1994, Shore 1998). These con 



ditions provide a rich testbed for the study of magnetic and 
hydrodynamical processes on these stars and conceivably on 
stars of related active O- and B-subclasses in which wind and 
photospheric line variability has so far gone undiagnosed. Sev- 
eral seminal papers on the He-strong star o Ori E have demon- 
strated that several peculiarities in the ultraviolet spectrum of 
this star can be explained by circumstellar clouds forced into 
co-rotation over the magnetic equator. For example, Groote 
& Hunger (1976) showed that high-level Balmer lines appear 
at certain rotational phases, demonstrating the presence of an 
occulting, comparatively low-density gas. This circumstellar 
medium can also be observed as periodically modulated Ha 



emissions (e.g.. Short & Bolton 1994) 



Using spectra from the International Ultraviolet Explorer 



(lUE) satellite. Shore and collaborators (Shore & Brown 1990 



(SB90), Shore et al. 1990) demonstrated that a common prop- 
erty of winds in both He-strong and He-weak Bp stars is that 
they are channeled by a dipolar magnetic field to a torus-shaped 
magnetosphere (cloud) located over the magnetic equator This 
torus may well be be located in a different volume than the 
plasmasphere suggested by radio studies (Havnes & Goertz 
Drake et al. 



1984 



1987 and Drake 1998), but it is likely to 



be related. Following the initial ideas of Shore (1987), Babel & 
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Montmerle (1997a) and Babel (1998) have discussed a model 
in which shocks are excited at the wind-torus interface, cre- 
ating high energy emission. The magnetic and rotational axes 
of these stars are separated by seemingly arbitrary angles. In 
most cases this angle is large enough to provide an opportu- 
nity to observe the time-variability of emission and absorption 
components of the strong lines and thus to study the azimuthal 
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structure of the clouds as the wind-cloud complex co-rotates. 
While these variations are often large and complex, all indica- 
tions are that they are rather stable over decades (e.g., SB90, 



Bolton 1994, Henrichs et al. 2001) 



Henrichs et al. ( |1998| ) recently proposed the existence of 
a distinct new class of magnetic rotators (hereinafter, rotating 
magnetic B stars). The three stars so far assigned to this group, 
PCep (HD 205021), HR6684 (HD 163472), and HD 184927, 
are characterized by redshifted emissions of the C IV resonance 
lines. To this small class, we would add the extensively stud- 
ied. He-strong star a Ori E. The Si IV and C IV resonance lines 
of this star also show redshifted emissions at phases of mag- 
netic extrema (Bohlender et al. 1986). We have examined the 



sequences of CIV and NV resonance profiles in three He- 
strong stars in the SB90 sample with known rotational peri- 
ods (HD 37017, HD 37776, HD 64740) and find that each ex- 
hibits phase-correlated redshifted emissions. An examination 
of A^£W5/P5-processed spectra of HD 36385 and HD 133518 
shows that the C IV resonance lines in these stars' spectra also 
exhibit variable, redshifted emissions. The periodic emissions 
of the NV and CIV lines of the hot 07 V star 9' OriC are 



well known. Walborn & Nichols ( 1994 ) and Gagne et al. ( 1997 ) 



have suggested that this star is a hot analog to the Bp stars. Be- 
cause of the similar variabilities of the red halves of their C IV, 
N V profiles, it is likely that these stars, and indeed perhaps all 
He-strong variables, should be included in this new class. 

In this paper we provide the first detailed analysis of spec- 
tral lines to determine conditions in co-rotating clouds attached 
to the rotating magnetic B stars. We undertake this study by 
using ultraviolet lines obtained entirely from archival echel- 
lograms obtained with the lUE satellite. We have three specific 
goals: (1) to estimate physical parameters of clouds (tempera- 
ture, density, extent, metallicity, and turbulence), (2) to study 
the relevance of Bp-like surface chemical composition (metal- 
depletion) on cloud properties, and (3) to examine the pro- 
cesses responsible for resonance line emissions and anomalous 
absorptions at certain phases. 

For these purposes we have chosen the first four stars dis- 
cussed above: o Ori E (the most well studied He-strong star), 
HD 184927 (a sharp-lined, magnetic He-strong star; Wade et al. 



1997)), pCep (a composition-normal star with magnetic wind 



properties similar to the He-strong stars; it is also the prototype 
of the p Cephei pulsational variables), and HR 6684 (a p Cephei 
variable with a moderate rotational velocity; Jerzykiewicz 



1972, Kubiak & Seggewiss 1984). These last two stars are pul 



sators and are not He-strong stars. The selection of these four 
stars was guided in part by our desire to include both intrinsi- 
cally slow (PCep) and moderately rapid rotators (oOriE and 
HR6684), as well as stars observed at both low and high as- 
pect angles. pCep and oOriE are probably intrinsically mod- 
erate and rapid rotators, respectively, viewed nearly equator- 
on. HD 184927 is an intrinsically slow rotator viewed at low 
inclination. Less is known abut HR 6684, but its low v sin / and 



possible rotation period near 3.75 days (Henrichs et al. 1998), 



magnetic and optical He I line variations (Wade et al. 1997). 
In order to describe the geometry of the magnetospheres, it is 
important that the sample of stars studied has a variety of ro- 
tational velocities and be viewed from a range of inclination 
angles. 

2. Reduction and analysis of archival lUE data 

2.1. Data reduction 



The lUE data for this program were the NEWSIPS-extracted, 
large-aperture high dispersion spectral files from the Short 
Wavelength Prime (SWF) camera. These were downloaded 
from the MAST [] data archive. The data used were all the 
lUE SWF/large-aperture images for oOriE (29), HD 184927 
(30), p Cephei (31), and HR6684 (41), except that for o Ori E 
the images SWF 04840 and SWF 27259 were excluded. The 
latter spectra were observed at isolated times and thus do not 
lend themselves to analysis of short-term variations. The high- 
dispersion SWF-camera data for both NEWSIFS and lUESIFS 
are known to depend on errors in echelle ripple correction and 
background subtraction. Ripple-correction errors were avoided 
by working with net flux data for individual ech elle or ders. To 
account for background extraction errors (Smith 1999), we de- 
rived corrections first in the p Cep and HD 1 84927 spectra by 
forcing the depths of the Si III, Si IV, CIV, and occasionally 
other resonance lines to zero. This procedure was generally in- 
sensitive to rotation because most of the resonance lines we 
studied reside well up on the flat or even damping parts of their 
curves of growth, and their line cores are well broadened by sat- 
uration. The zero-point corrections we derived were typically 
5-10% of the continuum level. For lines in short-wavelength 
orders, we interpolated corrections from the Ly a and other 
nearby orders. 



2.2. Computation of simulated spectra and spectral 
absorptions 

We used a suite of codes written by I. Hubeny and collaborators 
to compute model line profiles to compare with the lUE data. 
The first code is SYNSPEC, a spectral line synthesis code devel- 
oped for input no n-LTE or line-blanketed model atmospheres 
by Hubeny et al. (1994). We used standard LTE atmospheres 
by Kurucz ( 1993 ) in our work. SYNSPEC is embedded within 
an IDL wrapper (Hubeny 1996) such that we could easily cal- 
culate spectra with a variety of trial atmospheric abundances. 
SYNSPEC was then run to generate synthetic spectra similar to 
the observed ones by computing a disk-integrated, rotationally 
broadened spectrum and to convolve it to the spectral resolu- 
tion (13 000) of high dispersion lUE spectra. 

Our principal objective in this paper is to assess the ef- 
fect of torus-shaped magnetospheric clouds on a star's spectral 
lines (Shore & Brown 199C| ) as they alternatively move onto 
and off the projected stellar surface. To simulate the signatures 



suggest that it is a moderate rotator observed at low inclination, 
an interpretation supported by the single-wave variations of its 
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of a cloud on the composite spectrum, we used the Hubeny 
CIRCUS pvogTum (Hubeny |1996|, Hubeny & Heap |l996|). This 



code was written to compute line absorptions and/or emissions 
of a gas cloud situated either in front or off the limb of a refer- 
ence star. CIRCUS requires the user to specify physical cloud 
parameters such as temperature, density, geometry, composi- 
tion, microturbulence, column depth, and areal coverage factor. 
Although the code can accommodate as many as three sepa- 
rate clouds we considered only a single homogeneous cloud in 
nearly all our simulations. In its solution of the radiative trans- 
fer equation within a cloud, CIRCUS computes line emissions 
and absorptions separately. This feature can be used to compute 
the cloud spectrum either in LTE (adding differential emission 
and absorption terms) or in the two-level non-LTE approxima- 
tion, in which only the absorption terms are calculated. 



The code proceeds by first consulting a Kurucz (199C) 
line library of atomic absorption parameters and computing 
an opacity spectrum for a user-specified temperature and elec- 
tron density. The spectra were computed at a spacing of 0.01 
Angstroms, which oversamples the contributing lines. The op- 
tical depth in each line can then be determined using an in- 
put column density. The surface of the star is divided into a 
lOOx 100 grid and the local intensity spectrum is evaluated at 
each grid point, thereby taking into account the effects of fore- 
shortening and limb darkening. Typically, we computed CIR- 
CUS models in pairs. In the first case we assumed that a square 
cloud occults the stellar disk. In the second, we allowed the po- 
sition of the cloud to shift off the projected stellar disk so that it 
contributed only line emission to the spectrum. The dimensions 
of real clouds near a magnetic B star are actuall y like ly to be 
much larger than one stellar radius (Hunger et al. 1990, Bolton 



1994). Thus, we expect that our treatment underestimates the 



total effects of emission at non-occultation phases. However, 
since we modeled only differences in line strengths between 
occultation and non-occultation phases, this fact is not relevant 
to our analysis. 

The density and areal coverage factors of the cloud are nec- 
essary but not critical parameters in our analysis. For our initial 
models we started with a trial volume density of lO'^ cm^^, 
based on the appearance of high Balmer series member s dur- 
ing the occultati on pha ses of o Ori E (Groote & Hunger 1976 , 
Short & Bolton 1994). From the resulting column densities 
(§4-. 3) and an estimate of 67?* for the extent of the cloud in 



the line of sight from these authors, we derived a mean density 
of ^3-5 X lO'" cm^^*. As a compromise we chose a density of 
1x10^' cm^^, which we used in our models. The results of this 
study are not very sensitive to the electron density for the range 
considered. 

CIRCUS includes a provision for various doppler effects. 
The most important of these is stellar rotation, which CIRCUS 
treats by shifting the local intensity spectrum by a wavelength 
equivalent to the projected doppler shift at each point on the 
disk. The program computes the net doppler velocity between 
each projected element of the cloud and the background star 
along the observer's line of sight. In our model there are no dif- 
ferential velocities between regions of the cloud along the line 



of sight and the projected area on the star because the cloud 
is co-rotating. Additionally, because microturbulence can be 
an additional factor in determining the signatures of a cloud 
on a spectrum, we varied this as a free parameter. In general, 
values less than 20 kms^' resulted in similar cloud absorp- 
tion/emission spectra, so we used this value as a default. 

3. Basic picture 

To date only a few models of the dynamics of rotating magnetic 
B stars have been discussed for He-variable stars. Of these most 
work has focused on o Ori E. Borrowing from work of Groote 
& Hu nger ([l97^ [l982[ |l997| ; Hu nger & Groote [l993| ; Hunger 
et al. |1990 ), Shore and Br own ([l990| ). Babel & Montmerle 



( |1994 |1997^ , and Babel ( |1998[ ), this prototype and proba 



bly all He-variable stars exhibit variations because they have 
dipolar fields which are inclined to some extent with respect to 
their magnetic axes. There is now a broad consensus that closed 
magnetic loops from the two poles divert a weak wind emerg- 
ing from these surface regions. The particles emanating from 
opposite poles collide at high velocity and form a torus-shaped 
magnetosphere. The particles eventually leak out through the 
outer edge of the torus (defined by the Alfven radius) and also 
settle back to the star - to what relative extents is still unclear. 
The radiative wind theory of Springmann & Pauldrach 



(1992) suggests that there is a tendency in this region of T^^, 



log g space for H and He atoms to separate during their flow. 
This chemical fractionation is one of several complex feedback 
processes which determines the evolution of the wind. In our 
picture a relative preponderance of metal ions are directly ac- 
celerated by the magnetic polar wind along open field lines 
and leaves the star permanently. At the same time many helium 
atoms are decoupled early in the flow (as is hydrogen in cooler 
stars) and return back to the surface along roughly the lines 
of force that guided their upward flow. At the surface Bp-like 
chemical anomalies (He-richness, metal-deficiency) establish 
local vertical abundance gradients over time unless they are 
destroyed, e.g. through mixing. Wind particles guided along 
closed lines accumulate in the torus for a short time before 
they are lost at their edges. The length of time these particles 
spend in the torus is far shorter than the evolutionary lifetime of 
the star (e.g., Groote & Hunger |1982|, Havnes & Goertz Il984, 



Linsky et al. 1992), so the torus's composition will mirror the 



abundances found at the star's magnetic poles. 

The above description is meant only to lay out a broad pic- 
ture that is pertinent to this paper. Havnes & Goertz (1984) 
and Linsky et al. (1992) have discussed dynamic models of 



a magnetic plasmasphere in which nonthermal electrons emit 
gyro-synchrotron radiation in a current sheet some 10-20 R^ 



from the star (Phillips & Lestrade 1988). These models were 



conceived in part because of the abnormally high radio fluxes 
of these stars (Linsky et al. 1992| ). In contrast. X-ray observa- 
tions do not suggest an abnormal level of high energy emission 



(Drake et al. 1994 1. This same comment applies to the strength 



of the wind components of the ultraviolet resonance lines. 
Moreover, at a given rotational phase the variability patterns 
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oOriE 


HD 184927 


pCep 


HR6684 


Teff 


23 000' 


23 000^ 


27 OOO-* 


26000 


logg 


4.0 


4.0 


4.0 


4.0 


V sin i 


162-* 


152 


30 5 


100 


Prot 


1.1908106^ 


9.52961^ 


12.00092 ^ 


^^3.75** 


* = 0.0 


44241.537 


49987.388 


51238.15 


49953.03 



Table 1. Adopted Atmospheric Parameters. NOTES: (1) T^s 
is given in K, vsin/ in kms^', and the zeropoint of rotational 
phase (North magnetic pole crossing) as HJD + 2400000. (2) 
4> = 0.0 for HR 6684 might be quite inaccurate. (3) Sources are 
given as indices: 



1 Gro ote & Hunger |1997t 2 Wade et al. |1997[ 3 Gies & Lambert 
I99I 4 Bo lton et al. |1986| , 5 Campos & Smith |l980t 6 H esser 
et al. I1977L 7 Henrichs et al. bOOll 8 Henrichs et al |l998 



of the ultraviolet and optical lines appear very similar from one 
epoch to another at any given rotational phase. For these rea- 
sons the volumes emitting fluxes in these different wavelength 
regimes are probably not co-spatial. Thus, the physics of this 
emission is Ukely to be quite different. In this paper our philos- 
ophy will be to infer kinetic and thermodynamic properties of 
the magnetospheric torus-clouds within several stellar radii us- 
ing the ultraviolet and optical data alone. Additionally, we will 
be reluctant to invoke magnetic (or rotational) energy sources 
or sinks unless needed. 

Our examination of these questions proceeds by consid- 
ering first the obscurations of the photospheric spectrum by 
the torus-cloud. Since the clouds are locked into co-rotation 
with the surface, the photospheric and cloud contributions of 
a line profile cannot be separated according to differences in 
their shapes. Instead we must rely on the modulation of the 
cloud's absorptions when it occults the star once or twice each 
rotational period. Before proceeding to this analysis, we first 
discuss the relevant parameters required to compute the under- 
lying photospheric Une spectrum. 

Ultraviolet studies have shown that the He-strong stars are 
stars on or just evolving off the main sequence (e.g.. Shore 



199C) and hence have \ogg values in the range 3.5-4.0. We 
chose Teff values for our program stars from relevant papers in 
the literature and list them in Table |l]. Rotational periods and 
magnetic ephemerides are well known for these stars (except 
for HR 6684) and are listed as well. 



4. Spectral analysis of line absorptions 

4.1. The photosphere: metallicity and Teff 



Hunger et al. (1990) and Reiners et al. (2000) concluded that 
the metal abundances on the surfaces of o Ori E, are low at the 
magnetic poles and roughly solar-like over the magnetic equa- 
tor However, their analyses were performed on optical spec- 
tra with the implicit assumption that the effects of C II, Si II 
line absorptions from the magnetospheric clouds are negligible. 
The Reiners et al. metallic-line strength curves show the same 



offset in phase with respect to the magnetic pole crossing phase 
that our Fig.^ does, suggesting that enhanced optical and ul- 
traviolet absorptions are formed in the same locality. We tested 
this possibility by constructing additional CIRCUS models for 
the optical C II and Si II wavelengths and found that these were 
able to reproduce fractional equivalent width changes of 10- 
20%. Thus we are led to a contrary picture, namely that the 
variation of UV and optical lines and the continuum is due to 
absorptions within a co-rotating cloud. 

Our line synthesis analysis was conducted by first select- 
ing rotational phases for which line strengths in our lUE data 
were minimal and then using SYNSPEC to match spectra in the 
1853-1876 and 1906-1931 Angstrom regions. We chose pCep 
as a control star and discovered that the lines can be nearly sim- 
ilarly matched by model atmospheres having either a normal 
metal abundance and a high effective temperature (27 000 K) 
or a metallicity of one-tenth the solar value and a lower tem- 
perature (24 000 K). We were able to break this ambiguity by 
modehng the wings of the Si III A.1206 and Si IV A.1394, A.1403 
lines only with the higher temperature model. This discrimina- 
tion is possible because the ionization of Si^^ and hence the 
strength of A,1206 is more sensitive to Teg- than to metallicity. 
The high-temperature, normal-abundance model leads to less 
extended wings for X1206. Our adopted value Tgff = 27 000 K 
(and normal composition) agree s well with the value 26 600 K 
found by Gies & Lambert 1992 . Fig. ^ shows overplots of an 
lUE spectrum in the region of the Al III doublet with both the 
adopted and the 24 000 K synthesized models. 

Following the same procedure, we estimated the metal 
abundances for a minimum-absorption observation (SWP 
55488) of HR 6684 in the Al III and Fe III wavelength regions. 
As with p Cep, we found a trade-off between composition and 
stellar effective temperature, but this time both a low abun- 
dance and high Teff was needed to simulate the comparatively 
weak line strengths. As before, the weak wings of the Si III 
and Si IV resonance lines at the non-occultation phases con- 
firm this model. The fitting ambiguity was again settled in fa- 
vor of high temperature (26 000 K) but a low metal abundance. 
Fig.[l|b exhibits the comparison for both low ([Fe]= -1) and 
normal metallicty and this effective temperature as well as a 
low-metals, high temperature model. The metal lines for the 
normal-abundance, low-temperature case are of course even 
stronger For both effective temperatures the normal abundance 
model can be ruled out. For HR 6684 we estimate a a metallic- 
ity [Fe] = -1.0 ±0.3 dex (internal errors). 

We digress to point out that because helium-rich atmo- 
spheres have comparatively low continuous opacities, a synthe- 
sized spectrum matched to observed metal line strengths will 
tend to err on the high side of the true photospheric metallicity. 
However, such errors are likely to be small for the mild helium 
overabundances of 2-3 x as occur on most He-strong stars. A 
similar assessment can be made of the effects of magnetic fields 
to the measured line splitting, either due to Zeeman splitting or 
to the modification of the atmospheric structure. Such effects 
could produce in principle a line strengthening when the mag- 
netic poles are visible and thus would produce small-amplitude 
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Figure 1. a-d Comparison of lUE observations with ^KV^P/sC-synthesized spectra of the wavelength region surrounding the 
Al III doublet. The observations were made near the phase when the North magnetic pole crossed the stellar meridian: a for 
PCep, spectrum SWP 52432 is compared with SYNSPEC models having T^ff = 27 000 K, [Fe] = 0, and Tes = 24 000 K, [Fe] = 
-1, b for HR6684 (at line minimum; no field detection yet) spectrum SWP 55488 is plotted with models having Tgff = 26 000 K, 
[Fe] = 0, -1, and Tgff = 23 000 K, [Fe] = -1, c for aOriE, five spectra (SWP 07555-6, 15785-7) are compared with models having 
Teff = 23 000 K, [Fe] = 0, = -1 & -1.3, d for HD 184927 four spectra (SWP 14504, 1451 1, 55612, 55770) are compared with 
Teff = 26 000 K, [Fe] = -0.7 & -1.3; the latter model fits so well it is largely hidden by the observed spectrum. 



line strengthenings at this phase. Note that such effects go in 
the wrong direction to explain the general line-strengthening 
during the occultation phase. 

Our fits to lines of the He-strong stars aOriE and HD 
184927 were conducted by using the effective temperatures. 



23 000 K, from previous optical studies. For these stars the 
abundances refer explicitly to the magnetic polar caps. In 
Fig.[l|c we depict the synthesized spectrum in the AlIII line 
region with a low-absorption observation of oOriE at (j) w 
0.8, that is, between the crossings of the North magnetic pole 
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and the associated He-spot. For broad-lined spectra such as 
these, the accuracy of an abundance determination deteriorates. 
Nonetheless, our comparisons show a preference for the [Fe] = 
-1 ±0.4, or less. Thus, our results are in mild disagreement 



with the optical result with the Hunger et al. (199C) finding 



of a metal deficiency of only a factor of 2.5 but are compa- 



rable to the Reiners et al. (2000) results. The program star 



HD 184927, provides an optimal opportunity for line synthe- 
sis comparisons because its photospheric lines are virtually un- 
broadened. Fig.[l|d exhibits a match with an average of four 
lUE spectra in the AlIII line region with a metal abundance 
[Fe] w-1.3±0.3. 

We now consider line formations in the clouds, starting first 
with the weak metallic line absorptions formed at great dis- 
tances from the star and will then work our way inwards. 

4.2. Comparison of time -variable spectra 

Another approach to disentangling the spectra of the cloud and 
photosphere is by means of the Temporal Variance Spectrum 
{TVS). Its use opens the question of whether heterogeneous 
metal abundances across the star's surface contribute apprecia- 
bly to the observed variation of metallic line strengths. This 
spectrum permits a comparison of each line's activity with its 
mean strength. If photospheric abundances change from point 
to point on the star, the absorptions of all lines will be affected 
to some degree. On the other hand, if excess line absorptions 
occur in a cool circumstellar cloud at some distance from the 
star, low-excitation (and usually strong) lines are those we ex- 
pect primarily to show fluctuations in the TVS. 

Figs.^,e test this expectation for HD 184927. Because of 
its slow rotation this star offers the best case for detecting the 
variability of individual lines. In this figure we display two 
spectral orders containing prominent Al III and Fe III lines. Be- 
cause the excitation potentials of the lower levels of several 
other Fe III (and Ti II) lines are high, we give these values as 
well. The lower (TVS) plot shows that only those lines with 
excitations of a few eV or less exhibit significant TVS fluctua- 
tions. Fig.^ shows the same pattern for the AlIII line region 
of o Ori E. Clearly, all or nearly all the variations of the metal 
lines arise in a cool medium. Inspection of the line depths con- 
firms that they arise one-quarter of a cycle earlier or later than 
the transits of the star's meridian. 

To place limits on the variation of metals on the surface 
of HD 184927, we chose three lUE spectra obtained at faint- 
star (magnetic null) phases and compared line synthesis models 
containing weak, high excitation lines. Our best fitting models 
for the spectra in the two orders containing the Al III and Fe III 
lines suggested a metallicity of [Fe] = -1.0±0.3 (internal error), 
which is only marginally different from the metallicity of -1.3 
estimated from Fig. |l]d. The reader can discern the line activ- 
ity from the lower line in Fig.^, which is the TVS spectrum 
of all observations. Note that the Al III lines themselves still 
have a residual absorption, which we interpret as due to the 
cloud. Otherwise, the differences for most other (high-x) lines 
are generally nil. Plots for the spectrum in the Felll A-1914 re- 



gion show similarly small differences. We conclude that essen- 
tially all the UV line variations are due to absorptions from oc- 
culting clouds. Supporting evidence for this conclusion comes 
in addition from near-UV and visible-band data for o Ori E dis- 
cussed in the next section. 

As Fig.^ shows, the TVS profiles have double lobes. To 
investigate the cause of these structures, we averaged two sub- 
groups of spectra observed at phases of maximum and mini- 
mum strengths (hereafter referred to as clear-star and occulta- 
tion phases). The low-excitation lines in the clear-star spectra 
have triangular shapes and narrow cores in these plots whereas 
the occultation-phase spectra exhibit broad, rounded cores. We 
confirmed that these differences in shapes arise during particu- 
lar phases by inspecting individual spectra and seeing the pat- 
terns in them. Because the cloud-tori are the presumed products 
of a rapid deceleration of wind particles, it is natural to attribute 
the excess line broadening as a part of the deceleration process 
- a turbulence of 25-50 km s^ ^ Q 

4.3. The cloud spectrum and derived cloud parameters 

CIRCUS was used to compute the absorption of the stellar flux 
by the cloud (the so-called iron curtain) by subtracting the pho- 
tospheric spectrum, computed by SYNSPEC from the photo- 
spheric plus cloud spectrum computed by CIRCUS. This was 
done for a variety of trial temperatures, column densities, and 
occasionally other parameters such as volume density, micro- 
turbulence and areal coverage factor In general, the lower the 
cloud temperature, the higher the monochromatic absorption 
because of the increased influence of hydrogen opacity. Fig. ||a 
shows a comparison of the ratio of lUE spectra obtained at the 
maximum and minimum flux phases for a Ori E. The principal 
features of interest in this figure are the increased absorption 
below 1300 Angstroms and also near X1900. Our tests showed 
that these wavelength-dependent features can be reproduced by 
having cool temperatures, viz. Tiioud = 11500 - 13 000K. In 
Fig.|^ we show a fit with similar parameters to observations 
of p Cep taken from clear-star and occultation phases of the 
rotation cycle (but at similar pulsation phases). For both stars 
these cloud temperatures are consistent with radiative equilib- 
rium temperatures expected for material situated several stel- 



lar radii from the surface (e.g. Drew 1989 Millar & Marlbor 



ough 1999). We also experimented with two-zone cloud mod- 



els. Because inner-cloud temperatures are hkely to be higher 
in these models, and the opacities lower, the total cloud opac- 
ities at a given wavelength tend to be much as the same as for 
one-zone, low-temperature models. As a result, we found that 
for the weak Fe-curtain one cannot distinguish very well be- 
tween judiciously chosen multiple temperature and single low- 
temperature models. 



^ It is conceivable that this broadening could instead have a mag- 
netohydrodynamic origin, e.g. from the generation of Alfven waves 
by motions of field lines. However, for a dipolar field in a weak-field 
star like [3 Cep the velocities of such waves would probably be smaller 
than these values. Moreover, Alfven waves would produce motions 
transverse to the line of sight when the torus occults the star. 
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Figure 2. a-e Comparison of spectra in Al III doublet region for program stars; clear star and occultation pliases are shown in the 
top and lower spectrum (spectra offset for clarity). The TVS spectrum is shown at the bottom for the ensemble of spectra: a for 
P Cep, upper spectrum is from SWP 06234, 42939, & 52620, and the lower one is from SWP 46244, 46255, 52121, 52460, 52573, 
& 52594, resp., b for HR6684, upper spectrum is from SWP 14514, 55668, 55672, & 55732, c for aOriE upper spectrum is 
from SWP 07553, 07555, and the lower one from SWP 07534, 07536-7, 07560, 07583, & 07609, d for HD 184927 the upper 
spectrum (top) is from SWP 14468, 14475, 14487, 55669, 55685, 55701, and the lower one is from SWP 14504, 1451 1, 14681, 
55629, 55770, & 55781 (excitations are given for a few hues), e for HD 184927 as in d but for the wavelength range surrounding 
low-excitation Felll lines. 
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Figures, a and b: Observed and computed ratios of (quan- 
tity + 1) cloud absorption spectra for aOriE and pCep. a 
for a Ori E the observed curve (thick solid line) is the ratio of 
lUE spectra SWP 07555-6, 07558, 07560 = 0.40-0.59) to 
SWP 07583 ((|)=0.02, observed several hours earlier), binned 
one point per echelle order Dot-dashed and dashed line curves 
denote CIRCUS models with T,./„„,/ = 1 1 500 K and 13 000 K, 
with column densities of 6x10^^ and 1 x 10^^, and with cover- 
age factors of 100% and 50%, resp. The cloud spectra represent 
the fractional absorption plus one (in continuum units), b For 
P Cep the observations are taken from the ratio of SWP 52415, 
52594 to SWP 52514, 52653. The model has parameters T,./„„^ 
= 11 500 K, a column density of 1 x 10^^ particles cm^^, and a 
70% coverage factor 



As an additional test of the models for Fig.|^a, we reran the 
model s for the A-X3700-4000 range used by Groote & Hunger 
( 1976 ) to show that the torus occultation is accompanied by 
increased core absorption of the high-level Balmer-lines. Our 
models (not shown) exhibit similar line strengthenings as those 
observed (20-30%). This result further validates the interpre- 
tation that the line strength variations arise from absorptions in 
the cloud. 

To place our spectroscopic results in a context in which 
they can be compared with previous photometric studies, we 
summed the fluxes in four lUE echelle orders in the unin- 
terrupted wavelength range XA-1800-1905 and plotted them 
as a function of time. The resulting continuum (UVC) light 
curves for oOriE and pCep in Figs.^a,c (top curves) exhibit 
a pair of light dips corresponding to cloud occultations. The 
stronger of the two dips is 14% for oOriE and 11.5% for 
P Cep.^ We estimate the corresponding flux dips for HR 6684 
and HD 184927 to be 2±0.7% and 3^±1%, respectively, i.e., 
just barely detectable. Some interesting conclusions can be 
drawn about the chemical composition of the o Ori E cloud by 



The narrowness in phase of the dips for a Ori E may eventually 
turn out to be important in constraining the thickness of the torus in 
the direction of the magnetic axis, at least at several stellar radii from 
the star. 



comparing the amplitude of the primary dip in this star's UVC 
light curve, 14%, with the wl 1% variation in the Stromgren u- 
band (Hesser et al. 1977). To see what cloud conditions could 



make the lUE- and ground-based absorptions of o Ori E con- 
sistent with one another, we computed CIRCUS models in the 
A.X1800-1905 and XA. 3450-3650 (roughly, Stromgren-M) re- 
gions. Several cloud parameters were varied, including temper- 
ature, metallicity, turbulence, and column density. All our tri- 
als predicted larger absorptions for A-A-1800-1905 region than 
for the Stromgren-M region. The most important result of these 
tests was that the absorption ratio of these wavelengths was 
found to be highly sensitive to metallicity. For example, the ra- 
tio decreases markedly from 2.0 at normal metallicity to 1.4 
for [Fe] = -1. This sensitivity can be explained as follows: the 
fluxes of a B star in the M-bandpass are almost independent 
of metallicity because they are determined almost entirely by 
Balmer bound-free opacity. In contrast, in low-metal models 
the contributions of hydrogenic and metal line opacities are 
roughly equal in the A-A,1800-1900 region. Since the great ma- 
jority of the metal lines in this wavelength region are unsat- 
urated, their absorptions vary almost linearly with metallicity. 
Taking our computed ratio of 1 .4 for the low-metals case, these 
results suggest that observed M-band variations should be about 
1.4x0.11 = 14.5%, i.e. almost the same as the observed value 
of 14%. In contrast, a model with a normal cloud metallicity 
predicts a flux dip of about 22% for A,A-1800-1900, which may 
be excluded easily by the observations. Altogether, the metal- 
lic abundance of the cloud in a Ori E derived from a compar- 
ison of flux amplitudes in the middle- and near-UV is in very 



good agreement with our line synthesis analysis (§4.1). It is 



also in agreement with our general picture which requires that 
the chemical composition of the cloud should be essentially 
identical with the source material from the outflowing regions 
of the photosphere. 

4.4. Cloud parameters from weak resonance lines 

In this and the next subsection we consider the formation of 
the added absorptions of weak and/or strong resonance lines 
of silicon, carbon, and nitrogen by using CIRCUS models to 
match the added absorptions observed at occultation phases. 
We will also discuss the necessary roles of turbulence and non- 
LTE processes for strong resonance lines formed in the cloud. 

Table^ gives the extreme equivalent widths in Angstroms 
observed at clear-star and occultation phases. These strengths 
are measured with wavelength limits specified around line cen- 
ter (generally roughly given by ±.vsin i) in order to isolate spec- 
troscopic variations of a cloud from other circumstellar activ- 
ity, such as the wind. The technique of measuring line strength 
variations in specific wavelength bands of lUE spectra accu- 



rately was introduced by Shore & Adelman (1981) for the He- 



strong stars. To gain some perspective of how representative the 
variations are in the program stars, we also measured the equiv- 
alent widths of the same lines in the archival datasets of other 
He-strong stars studied by SB90 with known periods. We found 
that the range of variations for o Ori E is typical for this sub- 
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Figure 4. a-d Equivalent width indices of various lines of the program stars with rotational phase (periods given in text). In each 
case the index is defined as the reciprocal of the ratio of the fluxes within a AA- bandpass around line center (Table^ to the ratio 
of vahd fluxes outside these limits. The curves given are for Si IV A.1394, 1403 (mean), CIV A,1548, Si III X1206, CII X1298. 
The curves are renormalized by the factors shown in the annotations and separated by vertical shifts. 



class. For example, the range for the Si IV doublet of o Ori E is 
similar to that found for HD 37017, larger than for HD 37776, 
and smaUer than for HD 64760. 



For the weak resonance lines the differences between the 
observed equivalent widths in Table ^ in parentheses, is taken 
to be the total absorption contribution in the star's cloud(s). 
This amount can be compared directly with the detailed ab- 
sorption predictions of CIRCUS, given in Table I for pCep and 
HD 184927. For a stated variation of a given line, the star- to- 
star differences can arise from a variety of factors, but the dom- 
inant one is the star's metallicity. Thus, similar cloud absorp- 
tions for lines of the similarly metal-deficient stars HD 1 84927 
and HR 6684 will lead to similar derived cloud parameters. We 
also note that the dependence of the derived cloud parameters 
on metallicity is especially important for the strong resonance 
lines for which the greatest contrast exists between a (weak) 
photospheric line and a highly saturated cloud line. 



4.4. 1. C II, Si II and weak Si III lines 

In TableM we summarize our findings from Tables I and | 
by giving temperatures derived from CII, Sill, weak Si III, 
and NV absorptions. We kept microturbulence constant (20 
kms^', except for NV) in this compilation and took entries 
from the middle of the limits given in Table|. Before consol- 
idating these results, we discuss the status of the suitability of 
these lines in terms of possible practical problems. 

The CII A-A-l 335-6 doublet is a well-known temperature 
diagnostic for cloud-like conditions, but in lUE spectra the 
1335 A line is marred by the presence of an instrumental re- 
seau. In the spectra of oOriE the broadened wings of the 
1336A feature are blended by nearby low excitation CrII and 
Fe II photospheric lines, rendering A-1336 unusable for this star. 
Both lines are further marred by an interstellar component (par- 
ticularly for o Ori E). Despite these problems the C II lines no- 
ticeably broaden in occultation phase spectra of all our stars. 
Reference to TableH shows that the A,1336 variation is moder- 
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.545 (-) 


1.5 
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.417 
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0.8 
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.227 


.167 (.060) 


1.5 


.435 


.399 (.036) 


0.6 


.316 


.311 (.005) 


0.6 


.402 


.392 (.010) 


0.8 
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.459 


.339 (.120) 


1.5 


.479 


.416 (.063) 


0.6 


.524 


.467 (.057) 


0.6 


.406 


.394 (.012) 
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.353 


.198 (.155) 
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1.0 
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1855 


.551 


.322 (.229) 


1.6 


.633 


.319 (.314) 


0.8 


.772 


.449 (.323) 


1.3 


.629 


.353 (.276) 


1.4 
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.480 


.260 (.120) 


1.6 


.434 


.302 (.132) 


0.8 


.434 


.320 (.114) 


1.3 


.556 
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1.4 



Table 2. Observed Equivalent Width Extrema, NOTES: (1) Units of equivalent-width (EWs) and wavelengths are Angstroms 
were measured for AA, range noted. (2) Parenthesized entries after Min are differences between maximum and minimum EWs 
at midpoints of clear-star and occultation phases; a colon denotes an uncertain value. (3) Bracketed entries following the listed 
ion is the ionization potential in eV. (4) Blank entries for Si IV and N V indicate that models are not computed for these ions in 
Table || (5) N V entries (starred) for HD 1 84927 and P Cep refer to full-profile absorptions and are modeled in Table|^. 



ately large for HD 184927 and HR6684. Indeed, these are the 
only C II results that can be quoted with reasonable reliability. 

lUE/SWP spectra also include the Si II resonance lines at 
1264.7 A, 1265.0A, and 1309.3 A. The observed variations for 
the doublet are immeasurably small for all the program stars 
except oOri E. For p Cep, HD 184927, and HR6684 the small 
variations mean that the Si'+ population is low because the 
cloud temperature is high (or, less likely, very low) 

Variations in both weak resonance and low-excitation Si III 
lines provide convenient cloud temperature diagnostics. Con- 
sider first the forbidden Si III A,1892 resonance line, for which 
the log gf is 30000 times weaker than the A,1206 line's. This 
line is weakly visible in spectra of normal early B-type stars 
(Walborn et al. 1995), but it strengthens for our program stars 



during occultation phases. This is not surprising: in contrast to 
the photosphere, one expects a cool circumstellar gas with a 
long path length and a shallow temperature gradient to produce 
a strong absorption in a resonance line. 

The Si III Unes also shed light on the cloud density. Con- 
sider that the A,1892 resonance line in turn feeds the excited 
lower levels of the A.X1294-1303 multiplet of Si III. Both the 
observations (Table ^) and the CIRCUS models demonstrate 
that the variations of the Si III resonance and multiplet features 
are comparable and well correlated among themselves. These 



facts suggest that metastability of the ^P''3p (lower) level does 
not determine its atomic population and thus that the cloud den- 
sities are moderately high (> 10^' cm^^). 

To summarize this section, the ionization potentials of Si'+, 
C and Si ^+ (Si II, C II, Si III: Xip = 16-34 eV) together pro- 
vide an excellent bridge of diagnostics for conditions between 
the inner and outer cloud regions. From Table^ one can see 
that the line-formation temperatures decrease with ionization 
potential of the parent ion, as one would expect for a circum- 
stellar cloud having an outward-decreasing temperature. 

4.5. Resonance line absorptions as diagnostics 

4.5.1. Probing temperature and geometry 

According to our CIRCUS models, the ionization of sili- 
con and carbon most favors line formations at temperatures 
centered at 21 000 K and 26 000 K, respectively, for the ex- 
pected cloud densities. These temperatures are comparable and 
slightly higher, respectively, than what one might expect from 
radiative equilibrium models for plasma near stars having T^s's 
of 23 000-27 000 K. Thus although the resonance lines are not 
necessarily superionized features, they must be formed in a 
warm environment, thus close to the star As noted by various 
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Table 3. Computed Cloud Absorptions. NOTES: (1) Cloud temperatures are given in Kelvins, and microturbulences in kms^ . 
(2) Columns {CD 1-3) denote column density (units of 10^^ cm^^) for models. (3) Note that the absorption entries are computed 
from the same AA, limits (except for N V) as given in Table ^ (4) Pairs of entries close to observed values in Tableware italicized. 
(5) Computed EWs are often somewhat below or above the observed value. In such cases the table entries have trailing (-) or 
preceding (+) symbols, respectively. All models have coverage factors of 100%. 



authors (Fischel & Sparks 1980, Barker et al. 1982 Shore & 



Brown 1990, Henrichs et al. 1998), the strong resonance lines 



of CIV, Si IV and other ions undergo dramatic changes dur- 
ing the rotational cycle of each of our program stars. Fig. Q 
depicts the phase variations of an equivalent width index for 
several strong resonance lines for each program star.0 In each 



^ The depicted UVG curves for PGep and aOriE are time se- 
quences of binned fluxes constructed from the echelle orders extend- 
ing over A.A, 1800- 1905; errors are estimated to be ±3%. The PGep 
UVG data are large-aperture observations of 1991-4 and are corrected 



case this index is the ratio of the measured fluxes outside and 
inside the AA- wavelength limits for that line, as specified in Ta- 
ble ^ Similar indices computed from the blue and red halves 
of these wavelength windows give in all cases essentially the 
same fractional amplitudes. These checks indicate that the in- 
fluences of the wind at moderate velocities (or in the case of the 
C IV lines, the Fe III interstellar components) do not contribute 
appreciably to the variation of this index. The sequence of these 



for a 5% sinusoid found in the raw data at the pulsation period. The 
oOriE UVG fluxes are from 1980-1. 
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Line 


column 
density 


HD 184927 


pCep 


Cll 


1 ■ lO^-^ 


15 000 K 


(none) 


Sill 


3 ■ 10^2 


15 000 K 


16 500 K 


Si 111 


3 ■ 10^2 


16 500 K 


18 500 K 


NV 


1 ■ 1022 


30 500 K 


28 500 K 



Table 4. Synopsis of Best-Fit Cloud Temperatures for Lines of 
Ions Listed (from Table|[). NOTES: (l)Typical column densi- 
ties are in middle of best-fit parameter space. (2) Errors are typ- 
ically H-/-1000K except for C II(ISM and reseau flaws) and Si 
III (PCep). (3) Typical microturbulence values are 20 kms^', 
except for N V (50 kms^'). 



strong line variations may also be compared qualitatively to the 
variations of the Al III resonance doublet by consulting Fig. |l} 
The data in Fig.Q are well enough sampled to draw sev- 
eral conclusions concerning the geometry of the tori. First, the 
absorptions among different strong resonance lines peak at the 
same small ranges of phase. This fact suggests that these lines 
form over the same stellar longitude. However, for P Cep and 
a Ori E the ultraviolet quasi-continuum curves are offset from 
the central absorption phases of the strong resonance lines. This 
fact implies a war ping a t large distances from the star (see also 
Groote & Hunger 1997 GH97). Another conclusion from the 
sharpness of the maxima in these figures is that the heights of 
the clouds (perpendicular to their plane of symmetry) are at 
least one stellar diameter Otherwise, for short heights, their 
optically thick absorptions would be proportional to their pro- 
jected areas against the stellar disk. Such a geometry would 
impose an (unobserved) local minimum on the resonance-line 
absorption curve at the phase of mid-occultation. We used this 
fact to adopt an area coverage factor of 100% in our CIRCUS 
analysis of cloud parameters. 



4.5.2. Probing the wind 

To investigate the behavior of the blue-wing absorption from 
winds in these stars, we modified our equivalent width extrac- 
tion routine to compute flux ratios in a 2A segment of this 
wing to all valid fluxes in the order except those from the pho- 
tospheric line profile. With the contribution of the latter sup- 
pressed, we found no variation (i.e., to <10%) in the wings of 
either the Si IV or C IV lines with rotation phase among any of 
the program stars. Additionally, the blue edge velocities are sta- 
ble over the cycle and are also remarkably uniform from star to 
star. GH97 had already found that the coronal wind in a Ori E 
is phase-independent, and the constancy we find tends to ex- 
tend this conclusion to our other stars. This is a notable point 
because previous investigators have stated that the wind fluxes 
of HD 184927 vary from this star's magnetic pole to equator 



(Vauclair et al. 1991). This erroneous perception has been an 
important supporting tenet for the mo del of p hotospheric diffu- 
sive settling of helium (Michaud et al. 1987[ ). In addition, these 
observations suggest that the initial wind column, dubbed the 



jet by SB90, is not visible in these data, though indeed it must 
ultimately be visible at some level of photometric precision be- 
cause of the optical thinness in the far blue wings of the pho- 
tospheric lines. The difficulty of observing the jets could mean 
that their velocities are spread out over a large range. 

4.5.3. Probing density and turbulence 

Fig.|^ depicts the total ranges in absorption and emission from 
pairs of observations for the strong resonance lines for the pro- 
gram stars. Inspection of these figures shows that profiles at 
clear-star (dashed lines) phases differ from occultation phase 
profiles (solid lines) in the following ways: (1) the clear-star 
profiles of C IV and Si IV exhibit emissions in their cores and 
red wings (§||), (2) the Si III profiles exhibit narrower cores, and 
(3) the red and blue far-wings of these lines are strong relative 
to photospheric line models. 

As with the weaker metallic lines, one would hope to uti- 
lize the strengths of strong resonance lines at clear-star phases 
as baselines to compare with the absorptions at occultation 
phases. However, modeling the strong resonance lines is a 
problem complicated by the presence of noncoherent scatter- 
ing (partial redistribution) in the line source function as well 
as illumination of the torus by a diffuse radiation source (the 
jet). In general, it is difficult to separate the effects of the torus 
from both the photosphere and a jet. In particular, it is impossi- 
ble in practice to identify a reference photospheric profile from 
which to difference absorption and emission of the torus-cloud 
because the observed Si IV, C IV, and N V line profiles are pe- 
culiar at all phases, exhibiting strong absorptions during occul- 
tations and emissions at other times. Nonetheless, points 2 and 
3 above can provide hints of local cloud conditions. For exam- 
ple, point #2 suggests that the broad core of the occultation- 
phase profiles is a signature of velocity broadening. As to point 
#3, the extended, particularly red wings of the Si III and Si IV 
lines can best be understood as enhanced strengthenings from a 
dense cloud. Consider that our CIRCUS models for Si III A,1206 
and the Si IV resonance doublet indicate that densities less than 
10^' cm^^ cannot be distinguished from the clear-star (pho- 
tospheric) profiles. Our best models require cloud densities in 
the range of lO'^^'^ cm^^. Taking again the simplest paradigm 
that the cloud temperature increases inward toward the star, the 
wings of the strong silicon resonance lines suggest then that 
cloud densities increase inwards too. 

Traces of the variations so readily apparent in the Si IV, 
CIV, and NV lines are still present in Si III A-1206 and even 
in the Al III doublet. However, the difference between this line 
and the others is that the Si III line does not show any emission 
at any phase. While this circumstance might lead to the hope 
that one can model the A-1206 variations, we found in prac- 
tice that we could not reproduce the observed strengthening of 
this line without introducing a substantial microturbulence. As 
an example, consider the Si III line in pCep. Table ^ shows a 
cloud-added absorption of 343 mA. CIRCUS models s ugges t 
a maximum absorption for T,./w = 17 000-19 000 K (§ M.l| ). 



Even for column densities of 3-10 x 10^2 cm ^ and a turbu 
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Figures, a-d Comparison of spectra of Silll, SilV, CIV, and N V resonance lines of minimum (clear-star phase; dashed line) 
and maximum (occultation phase, solid line) absorption strength (shifted for clarity). Except for p Cep the normalized TVS 
spectrum is also given, a for p Cep, the coiTesponding spectra are from SWP 57573 (dashed line) and SWP 52514 (solid line), b 
for HR 6684 the corresponding spectra are from SWP 55732 and SWP 55526, c for o Ori E the corresponding spectra are from 
SWP 07560 and occultation phase SWP 15787, d corresponding spectra for HD 184927 are from SWP 55652 and SWP 55685. 



lence of 20 km s^\ our simulations predict an absorption in the 
range 138-259 mA. Low-temperature, pure-absorption {non- 
LTE) models also cannot reproduce the observed strengthen- 
ing. If one introduces a high turbulence of 50 kms^' and also 
adopts a reasonable temperature and column density (17 000 K 
and 3 x 10^^ cm^^), the observed absorptions («340 mA), can 
then be reproduced. The Si III line's sensitivity to turbulence 
implied here is a consequence of the saturation in the core. 
One may assume that the broad cores of the CIV and Si IV 
lines at occultation phases are also due to turbulence. As for 
Si III A-1206, moderately high turbulences and moderate tem- 
peratures produce good fits to the AlIII doublet absorptions. 



noted in Table |[ We also note similar good fits for HD 184927. 
This is not surprising in view of the similar ionization potential 
of this ion: the Al III lines are probably formed in nearly the 
same region of the cloud as Si III A,1206. 

4.5.4. The N V absorptions 

The anomalously strong absorptions of the N V 1238A, 1242A 
doublet at occultation phases first drew our attention to the ro- 
tating mag netic B stars. According to stellar atlases (e.g. Wal- 
born et al. 1995), the spectra of normal B0-B2 V stars show 
no trace of this doublet. Even for the hot star 10 Lac (09; Teff = 
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35 400 K), we estimate that this line has a photospheric strength 
of only 175 mA (see the Brandt et al. |1998| atlas). Our SYN- 



SPEC models reveal that their presence is difficult to discern in 
photospheric spectra in stars with Tiff's below about 28 000 K 



because of nearby line blends. Grigsby & Morrison (1995 ) have 



recently noted the anomalous strengths of N V lines in spectra 
of the so-called nitrogen enhanced late-O stars (Walborn et al. 



1985). They attribute them to a photoionization of the wind by 



X-rays. However, the wind explanation cannot apply to the ro- 
tating magnetic stars in our sample because the lines do not 
show extended, high-velocity (blue) wings. Moreover, the N V 
strengths in our program stars tend to be larger than the features 
in these O-stars. By elimination, we conclude that the N V ab- 
sorptions are formed in the circumstellar cloud(s). 

The N V absorptions in these stars are so large («750 mA 
and «530 mA over thtfull profile for P Cep and HD 184927, 
respectively) that they cannot be caused by enhanced abun- 
dances or strong winds, or fit by our conventional strategies. 
For example, our standard cloud models, with both high cloud 
temperatures (~33 000 K) and a very large high column density 
(10^^ cm^^) can not reproduce them. It is only when one con- 
siders models with both high cloud temperatures and microtur- 
bulences together that the strong absorptions can be matched. 
For such combinations the precise strength of the absorption 
is insensitive to column length because the line is already so 
opaque. To assess the effects of increasing both these parame- 
ters, we carefully computed the underlying photospheric spec- 
trum, taking into account the star's Tgff, \ogg, and vsin/. The 
CIRCUS models showed that the N V line is extremely sensi- 
tive to temperature, a fact that greatly restricts the range of vari- 
ables in successful models. For P Cep we settled on a promis- 
ing pair of models (normal abundances, column density of 10^^ 
cm^^, microturbulence = 50 kms^'). Then, a change from 
28 000 K to 29 000 K strengthens the equivalent width from 
585 to 742 mA. Increasing the microturbulence from 50 to 100 
kms^' produces exactly the same increase. Thus, both micro- 
turbulence and temperature are sensitively constrained by the 
strengthening of this line. We obtained similar temperatures 
for oOriE and HR6684, with a turbulence of 50 kms^'. For 
HD 184927 we were able to match the absorption with a tem- 
perature of 30 000-3 1 000 K, assuming a column density of 1- 
2x10^^^ cm^^ and 50 kms^'. The best fit models are summa- 
rized in Table 

Fig. H exhibits a model fit of the N V doublet for p Cep at 
the occultation phase. The fitted cloud absorption profile was 
computed with a CIRCUS model having a coverage factor of 
100%, heated T,./„„,/ = 29 000 K, a column density of 3.2 x 10^^ 
cm^^, a microturbulence of 50 kms^', and a further con- 
volution with a 100 kms^' gaussian. This macroturbulence- 
smearing was necessary to remove the corners of a boxy-shape 
profile produced by microturbulence alone. The parameters re- 
quired to approximate the N V absorption of HD 184927 are 
similar Because HD 1 84927 has a low surface metallicity, the 
required cloud temperature is a little higher, Ri3 1 000 K, and 
the derived micro- and macro-turbulences are then 50 kms^' 
and 100 kms^', respectively. 
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Figure 6. Comparison of mean spectrum of the N V doublet 
for P Cep compared to a CIRCUS model fit (upper panel). The 
observations are taken from the the occultation-phase spectra 
listed in Fig.^. The fitted model represents a cloud with the 
temperature and microturbulence shown. The dashed fit is the 
result of smoothing with a a 100 km s^' gaussian. Lower panel 
shows the TVS spectra divided by mean spectra for the ensem- 
ble of large- (LGAP) and small-aperture (SMAP) observations 
of P Cep. 



5. Resonance line emissions 

The spectroscopic hallmark of the rotating magnetic B stars is 
the strong variation of their C IV and N V doublet lines with 
rotational phase, ranging from strong absorption to redshifted 
emission. Since these variations are not necessarily well or- 
dered, only a few generalizations can be made about their be- 
havior One overriding general characteristic is that the emis- 
sion and/or absorption components are always present to some 
degree. A second characteristic is that although small differ- 
ences can arise in same-phase profiles of different epochs, 
the general variations are remarkably repeatable over at least 
several years. For example, a narrow redshifted emission fea- 
ture seems to be a common feature in the 1238 A profiles of 
HD 184927 at phases near 0.85 and 0.15, i.e. at the begin- 
ning and end of the range of clear-star phases. For some stars, 
such as HR 6684, the N V profiles almost always show emis- 
sion components, even when much of the profile is strongly 
in absorption. For other stars like a Ori E, the N V emissions 
are difficult to detect at all. The weakness of N V emissions in 
O Ori E, along with the weakness of Si IV emissions for all of 
the program stars, probably explains why these emissions have 
gone unreported for so long in studies of He-strong stars. 

All four of our stars show detectable emissions in lines of 
at least two ions. In two of our four program stars, the N V lines 
exhibit the strongest emissions, while for HR 6684 or a Ori E it 
is the C IV emissions which are the strongest. 

In some cases, the N V emissions of these stars are almost 
completely confined to the red side of the profile, and indeed in 
all cases the TVS is larger there. The lower portions of Figs.|^ 
b, c, d and ^ depict the TVS spectrum divided by the mean pro- 
file for the C IV profiles or N V profiles. We depict the ratio of 
the TVS to the mean spectrum in order to account for absorp- 



Myron A. Smith & Detlef Groote: Wind circulation in selected rotating magnetic early-B stars 



15 



tions of the strong line core which would otherwise suppress 
fluctuations at low velocities. The TVS-ratio profiles generally 
peak at +60-125 kms^', and their red wings extend to +300- 
400 kms^' in each star with remarkable uniformity. Thus, the 
edge velocities are not related to the peak emission strength, 
the star's metallicity, pulsational status, or vsini. 

A comparison of the TVS-ratio spectrum for P Cep among 
both large- and small-aperture observations in Fig.|| gives a 
qualitative idea of the robustness and constancy of the activ- 
ity spectrum over even widely separated epochs for this star. 
It is interesting that the TVS profiles for C IV and N V do not 
exhibit a tail extending to the far blue wing. This fact shows 
that the stars' winds do not vary much over stellar longitude. 
In all, one finds that the redshifted emissions are uncorrected 
with wind visibility in the blue wings of the lines. Likewise, 
the emissions do not appear to arise from the wind, so one must 
look to a different mechanism to explain them too. 

We can quantify the temperature range required to repro- 
duce the N V emissions by again using CIRCUS models, again 
taking p Cep as an example. As one increases the cloud plasma 
temperature in these models, the N V lines first show increased 
absorptions at temperatures. These absorptions peak in the 
range -35 000-40 000 K, both for LTE and non-LTE models. 
As the temperature increases the line goes into emission. One 
may place interesting upper limits on the line formation re- 
gion of the N V emissions because at temperatures in excess 
of 40 000 K these regions of the cloud would contribute an un- 
acceptably large fraction to the combined continuum flux of 
the star and cloud. Such models would produce a nearly sin- 
uosoidal modulation in continuous flux as the torus co-rotates 
around the star. This is not indicated by our continuum light 
curve data, so we can conclude that any continuum flux from a 
bright emission patch near the poles is negligible. An additional 
consequence of these arguments is that the N V line emissions 
are optically thin. This result is also supported by the high ra- 
tio (1.5-2) of the relative r.m.s. amplitudes of the two lines in 
Fig. H (the doublet components have a gf ratio of 2). 

In assuming the optical thinness condition and an areal cov- 
erage factor of unity, it becomes straightforward to model the 
A,1238 emission equivalent width. Again taking p Cep as an ex- 
ample, the observed emission over the full profile is « -220 
mA. In the models there is only a small trade-off between 
gas temperature and column density, e.g., a choice of param- 
eters 49 000 K and 3xlO" cm^^ or 50 000 K and lO'** cm^^; 
the limiting condition for X;^i238 = 1- Because of the sensitivity 
of N V ionization with temperature, the results for HD 184927 
(full profile equivalent width — -0.055 A) are similar to those 
for p Cep: we find an optimal temperature of 45 000 K for the 
resulting length of the NV emission region is 1x10^' cm^^, 
the N V-emitting region comes out to be prohibitively thin, only 
— 100 km. This is roughly lOOx smaller than a scale height in 
the torus. The shock depth can be reconciled to the expected 
range of values if the assumed coverage factor for N V emis- 
sion is <1% rather than the initially assumed value of 100%. 



6. Interpretation 

A key result of the previous section is that the UV line ab- 
sorptions of a rotating magnetic B star vary because of a cloud 
passing in front of the star Thus, there is no need to invoke (nor 
evidence for) an inhomogeneous metallic composition on the 
surfaces of these stars. The observations also suggest that the 
mass in these clouds is nearly constant over time. To maintain 
a constant mass, the wind accreted by the cloud must somehow 
settle back to the star and/or escape through its outer point. In- 
deed, perhaps it is the settling of cloud material onto the surface 
that is responsible for the high densities we found in the inner 
cloud regions. As a consequence of this recirculation, the cloud 
has virtually the same metallic composition as the star's polar 
caps and indeed the rest of its surface. 

The numerical analysis of line strength variations in §0 in- 
dicates that the cloud line components form in an extended en- 
vironment of a column density of roughly 1 x 10^^ cm^^. In the 
detailed models of the cloud spectra of HD 184927 and p Cep, 
the lines are formed over a range of temperatures consistent 
with the stellar photosphere's radiation field. Although there is 
no evidence for long-term changes in cloud conditions, there 
is a strong suggestion of small- and large-scale turbulence in 
the inner regions. The strong wings of the Si IV and Si III res- 
onance lines argue for increasing cloud densities in these re- 
gions, which seem to approach photospheric values at the inner 
edge of the torus. 

The redshifted emissions in the NV and CIV lines ob- 
served during the magnetic pole-on phases demonstrate that 
downflows occur in the vicinity of the rotating, magnetic B- 
stars. In our picture they are most likely to arise from shocks 
because of the rapid deceleration of the wind particles as they 



impact the pre-existing torus (cf. Babel & Montmerle 1996) 



This produces a redshift during the phase when the magnetic 
pole is visible because the wind has followed the curved lines 
of force and is then moving away from the observer However, 
additional possibilities are possible. One should also consider 
the possibility of a stalled wind which impacts the star's sur- 
face near its point of origin. (In principle, it is also possible 
that magnetohydrodynamic al proc esses, or Kelvin-Helmholtz 
instabilities; see Shore et al. ( 1987 ), operate, but we see no need 
for this in view of of the potential of steady-state hydrodynam- 
ical iteractions alone.) 

In examining the evidence, we have found two emission 
characteristics that cast doubt on whether the emissions are 
produced by reaccretion, and we review them in the follow- 
ing. We discovered the first of these characteristics by recall- 
ing that SB90 had drawn attention to red-shaded C IV emis- 
sions in two He-strong stars for which periods are not known, 
HD 96446 and HD 58260. SB90 also noted that these features 
cannot arise from spherically symmetric shells, i.e. they are 
not true P-Cygni profiles. The SB90 star HD 96446 is likely to 
meet our criteria for membership among the rotating magnetic 
B stars because it also shows NV line emission. In addition, 
its emissions are variable. The case of HD 58260 is particu- 
larly interesting in this analysis for two reasons. First, SB90 
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have pointed out that we probably observe this star at a pole- 
on aspect (the magnetic and rotational poles are nearly coinci- 
dent), so its emissions do not vary. Second, this star's C IV line 
emissions are both optically thick and strong - indeed, it at- 
tains a level of 0.51t„,„. In attempting to model C IV emissions, 
we found that it was not possible to reproduce this emission 
strength within the context of our models. Even our best model 
parameters (T^-ii^ni = 30000 K, turbulence = 50-100 kms"', 
high column densities, coverage factor = 1) produce C IV emis- 
sions of only 0.20 leant- We also found that the efficiency for 
producing the emission peaks very sharply at 30 000 K, so areal 
coverage factors of at least 3-4 x are probably necessary. Since 
this factor is significantly greater than one, the emissions can- 
not be produced just on the stellar surface. The area of the 
cloud-torus is more than adequate to produce the total emis- 
sion, so we are led to the conclusion that the emissions are 
caused at the wind-cloud interface. 

The second characteristic we used to decide between the 
competing shock models for line emission came from the de- 
tailed examination of the evolution of the C IV and N V emis- 
sions in both p Cep and a Ori E. We discovered that an emission 
component appears first at the end of the occultation phase. In 
a model in which the downflows fall onto the star's surface, 
we would expect the line-emitting material to share the star's 
rotation because of its co-rotation with the star. As the cloud 
clears the receding limb of the star at the end of its occulta- 
tion, emission is produced as the failed-wind stream shocks 
at a magnetic pole (now appearing at the opposite limb of the 
star). The emissions from this region will have an algebraically 
smaller radial velocity both because of the star's rotation and 
because the projected component of the infall velocity at this 
time is nearly zero. The emissions should then migrate red- 
ward smoothly across the line profile in this model as both 
velocities become more positive. The actual observations in- 
dicate that the emission component first appears weakly at a 
common point in the line profiles of both stars (^^-200 km s^') 
and then migrates to the red. This smooth behavior is suddenly 
interrupted by the strong appearance of a red emission compo- 
nent at H-O-lOO kms^'. This particular evolution is difficult to 
understand as a consequence of shocks occurring on the stel- 
lar surface, particularly for P Cep and a Ori E, which have very 
different rotational velocities. In contrast, the disk interpreta- 
tion is almost too easy, suffering only from an embarrassment 
of several nonunique geometrical interpretations. 

A related question is how the N V absorptions are formed 
during occultation phases, sometimes even coexisting with the 
emissions. The hypothesis that the NV emissions are emit- 
ted from the torus interface provides an advantage of restrict- 
ing the formation of the emission to the same general site 
as the absorption. Both sites require a high temperature. Un- 
like the emission component, the absorption is optically thick. 
CIRCUS models show that absorption of N V is necessarily 
formed in a turbulent medium having a warm temperature of 
about 30 000 K. These results suggest that the N V absorption 
is formed in a large region of the cloud in which C IV and N V 



are overionized. We suggest that the wind-cloud shocks gener- 
ate this overionization. 

Still, it is not clear whether the heating required to produce 
the N V absorption occurs from a conversion of mechanical or 
radiative energy. The problem with either mechanism is that it 
is difficult to heat a large enough cloud volume to produce the 
amount of absorption observed. For example, turbulence can 
circulate shock-heated gas but probably only to a few scale 
heights. For representative cloud parameters (30 000 K, lO'^ 
cm^^) this length is probably only 10^ km or less. Radiative 
models encounter similar difficulties if the radiation from the 
shocks is emitted primarily in the EUV The photoionization 
radius for in our cloud models is only '-^lO^ km. In con- 
trast, CIRCUS models that match the N V absorptions require 
a full areal coverage factor, so the scale length for energy dissi- 
pation must be at least IR*. The radiative model can be saved 
if the shocks generate soft X-rays (1/2-1 keV) which would 
be transparent to the cloud over its full length. Such irradiation 
would produce a broad distribution of ion states (Cohen et al. 



1994 ) and would thus account for the observed strong absorp- 



tions produced by the cloud over a range of ion stages. The ev- 
idence for X-ray heating on these stars is ambiguous. Although 
the mean X-ray levels are known to be low, little is known yet 
about variability. The kinetic model in turn might be reconciled 
with models of relativistic particles penetrating to great depths 
into the torus of aOriE. Alternatively, it may be that porosity, 
rarefied zones, and/or Alfven waves play roles in transporting 
mechanical energy into the cloud. Although we favor the radia- 
tive model, all of these possibilities are speculative. Thus, the 
question of the origin of the N V absorptions remains open. 

7. Conclusions 



The present work confirms the finding by Henri chs et al. ( 1998| ) 
that a separate class of rotating magnetic stars exists among the 
early-type main sequence B stars. Taking pCep as a case in 
point, one finds that a dipolar field of only a few hundred gauss 
can suffice to guide a polar wind from the surface to its equato- 
rial regions via closed magnetic loops. The wind particles are 
abruptly halted near the magnetic equator where they establish 
a quasi-static, torus-shape structure (SB90). Our picture differs 
from the picture of GH97 and other investigators in two im- 
portant respects. First, we have found that variations of the UV 
and probably optical lines (other than of helium) are due to ab- 
sorptions in the torus-cloud. Second, we have remarked that the 
blue wings of the Si III, Si IV, and C IV resonance lines remain 
remarkably constant with rotational phase for all the program 
stars. Let us examine a few ramifications of these results. 

If indeed the variable line absorptions are caused by the 
clouds, and if the metallic composition is nearly constant across 
their surfaces, their low-metallicities all over the surface will be 
maintained only if cloud material settles back to the star more 
efficiently than its wind streams upward. One may speculate as 
to how this might happen. In our model, the result of this set- 
tling is that the inner torus has a high density and thus controls 
the configuration of the local field lines, rather than the other 
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way round. For example, the weight of the particles over the 
equator may cause the field lines to sag, as in some dense so- 
lar quiescent prominences, and particles could return directly 
to the surface. Turbulence may also help recirculate cloud par- 
ticles to the star. 

The phase-independent nature of the wind is more diffi- 
cult to understand. GH97 pointed to the two types of wind, 
one originating from the surface (magnetic pole), the second, 
an exo-magnetospheric wind originating at the outer edge of 
the torus. If the absorption contributions of the two wind com- 
ponents were equal in the blue wings, one would expect the 
fluxes of their blue wings to vary sinusoidally twice per cycle - 
yet this is not observed. However, even if the wind mass fluxes 
were comparable, it is likely that their observable contributions 
in the blue wings of the line profiles could still be unequal for 
a few reasons. For example, the exo-magnetospheric wind will 
be more likely to remain visible over a wide range of phases, 
both as a consequence of its inverse-square density dependence 
with distance and because its streamlines will arc backwards to 
maintain angular momentum conservation. These curved tra- 
jectories will cross our lines of sight over a variety of rotational 
phases as they go off to infinity, diluting any modulation of 
wind absorption at any particular phase. The key requirement 
of all such arguments is that the wind remains visible over long 
path lengths. 

As outlined in §^ a necessary requisite of our picture of 
rotating magnetic B stars is that closed field loops divert the 
wind flow from the magnetic poles to the torus at the mag- 
netic equator. In this picture stalled wind particles return to the 
magnetic caps and establish a large composition gradient. This 
does not occur in nonmagnetic stars because envelope mixing 
destroys composition gradients before they can be established. 
Also, similar processes may take place in hot Bp analogs (e.g., 
6' OriC), but its wind is so strong that chemical separation 
cannot be expected to take place. The results of this paper sug- 
gest indirectly that composition gradients are destroyed in mag- 
netic, pulsating stars, such as in p Cep and HR 6684. Yet, there 
are fundamental puzzles in even this simple picture. For ex- 
ample, what is to be concluded about HR 6684, which has low- 
metal but normal-He abundances? Perhaps the metals leave this 
star in the wind while both hydrogen and helium atoms are 
reaccret ed. Su ch circumstances are envisioned by Hunger & 
Groote ( 1999 ), but only in late-type B stars. A second puzzle 
is the Be-like episodes e xhibit ed by pCep perha ps eve ry 30 
years (Kaper & Mathias |l995i Panko & Tarasov |1997| ). The 
tight scatter in the resona nce lin e absorptions in Fig.|4|a and in 
Fig. 5 of Henrichs et al. ( ZOOl ) suggests that the presence or 
absence of a flattened Be disk does not affect torus-cloud sig- 
natures noticeably. Thus, if these ejections are related to the 
Be-phenomenon, the UV phenomenology of the two types of 
circumstellar structures is vastly different (in addition to their 
physical morphologies). 

In §1 and §5 we listed 13 probable members of the rotat- 
ing magnetic B-star group (including 9' OriC). Rotating mag- 
netic B stars appear also to include an indeterminate number 
of early-B main sequence stars with conditions that promote 



envelope mixing (pulsation, rapid rotation). One such example 
may be the unusual star HD 144941 (see also GH97). This star 
lies within the region of the H-R Diagram populated by the He- 
strong stars and shows extreme He-richness and metal deficien- 
cies. Its high galactic latitude, and high radial velocity (-52.6 
kms^') led Harrison & Jeffery (1997a 1997bi and others to 



consider it as a member of the hydrogen-deficient B star class. 
Moreover, Harrison & Jeffery ( 1997b| ) found that the Fe abun- 
dance is just as deficient as the CNO elements and suggested 
that these abundances are primordial. To test this result, we in- 
spected two high-dispersion lUE spectra obtained for this star 
reveal redshifted emissions in both the C IV and Si IV lines as 
well as a strong wind. These spectra show anomalously strong 
N V absorption so perhaps HD 144941 is a member too. We 
must also remark that the search for magnetic and rotational pe- 
riods is very labor and telescope intensive, so new members of 
this class might be found in the long run by examining spectra 
of early-B stars for the presence of anomalous N V absorptions 
and/or Si IV, C IV, N V emissions. 

The physical attributes of the tori around the rotating mag- 
netic stars can also be used to clarify our understanding of clas- 



sical stars. Several studies (see Balona ZOOO) have suggested 



that co-rotating clouds are responsible for periodic spectro- 
scopic and/or photometric varations in these stars. The diffi- 
culty in distinguishing between this explanation and nonradial 
pulsation is that optically thick clouds can mimic cool regions 
on the atmosphere. We are optimistic that some of the results in 
this paper that indicate the presence of redshifted components 
in resonance lines, the selective variation of low-excitation 
lines, and high turbulence together can help to discriminate be- 
tween models responsible for the variability of these stars. 

A general understanding of physical conditions leading to 
co-rotating magnetospheres will be aided by future studies of 
the magnetic rotating B variables (with the probable extension 

[T990| ), hot He- 



to the sn He-weak stars (see Shore et al. 1987 



strong analogues like 0' OriC, Ap stars like IQ Aur (Babel & 



Montmerle 1997b), magnetic Herbig Ae-Be stars, and late-type 



active stars such as AB Dor The results of this paper point to a 
number of desiderata which will clarify the nature of the class 
of rotating magnetic B stars in particular and also contribute to 
its global understanding. We close by listing a few of them; 

- the lUE archive can be conveniently investigated to search 
for additional examples of early-B type stars with anoma- 
lous N V absorptions and redshifted emissions among its 
C IV and N V resonance lines, 

- far-UV snapshot observations can be undertaken to deter- 
mine whether the anomalous absorptions and redshifted 
emissions extend to the O VI doublet, 

- the X-ray fluxes from these stars can be monitored around 
their rotational cycles to shed some light on the production 
mechanism of the C IV, N V absorptions, 

- optical observations of the high Balmer lines can be carried 
out over rotational phase among all the candidate members 
of this class with a view of establishing whether cloud den- 
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sity varies with time and as a function of relevant stellar 
parameters, 

- spectral imaging can profitably be carried out in the Ha Une 
and ultimately in the strong ultraviolet resonance hnes from 
space-borne interferometers. 
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